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مقایسه‌ی کاربرد مدل‌های آشفتگی 511 و 1.۳5 در پیچ تند * 


جواد مظفری(٩‏ امیر صمدی ۲ سید اسداله محسنی موحد(" .. داود داودمقامی* 


چکیده شناخحت الگوی جریان و مکان‌های ف رسایش و ر سو بگذاری در پی چآبراهه‌ها دارای اهمیت زیادی می‌باشد. اما استفاده از مدل‌های 
فیزیکی و برر س یآزمایشگاهی الگوی جریان در پیچ‌ها دارای محدودیت است و مدل‌های ریا ضی می‌توانند کمک شایئی به شناحت الگوی 
جریان در پیچ‌ها کنند. در این تحفیق مدل‌ها یآشفتکی 8511 و 5ظا در پیچ تند ی کآبراهه‌ یآزمایش‌گاهی با استفاده از نرم‌افزار فلوثنت 
مقاي سه شده‌اند. آبراهه‌ی آزمای شکاهی, کانالی با پیج تند | ست که د رآزماٍ شگاه هیدرولیک داد شگاه ,31 م ستقر می‌با شد. برای برر سی 
مدل‌ها یآ شفتکی, از پارامترهای متو سطگیری شده در عمق | ستفاده گردید. هم‌چنین برای برر س یکمی مدل‌ها یآ شفتکی» تعدادی پروفیل 
سرعت طولی اتتخحاب شد و با پروفیل‌های اندازه‌گیری شده مقایسه گردید. نتایج بررسی سرعت طولی متوس طکیری شده در عمق نشان داد که 
مدل 15 توانایی تعیین دقیق نقاط مهم جریان را ندارد. این مدل, منطقه‌ی جداش دگی جریان را در زاویه‌ی ۸۵ درجه و با بیشینه عرض ۶۰ 
درصد نشان می‌دهد در حال ی‌که در اندازهگیری‌ها, این منطقه دارای بیشینه عرض *1 درصد عر ض کل و در حدود زاویه‌ی ۷۵ درجه واقع 
می‌باشد. در مقابل مدل 5 به صورت منا سبی منطقه جداش دگی را ذشان می‌دهد و موقعیت‌های بيشینه سرعت طولی مت و س طگیری شده در 
عمق را نیز منا سب‌تر از مدل 15/۷ م شخحص می‌کند. برر سی سرعت عر ضی متو سطگیری شله در عمق نیز دشان داد که موقعیت‌های تغییر 
جریان عرضی به‌وسیله‌ی مدل 185 مناسب‌تر از مدل 85۷ شناعته شده است, نهایتاً مفایسه‌ی کمی پروفیل‌های سرعت نشان داد که مدل 1۳85 
حطای کمتری دارد و نتایج مناسب‌تری را ارائه می‌دهد. 


واژه‌های کلیدی پیچ تند. فلوئنت. سرعت متوسط گیری شده در عمق مدل‌آشفتگی. 
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تاریخ دریافت مقاله ٩۳/۵/۲۱‏ و تاریخ پذیرش آن ۹۶/۳/۲۲ می‌باشد. 

(۱) استادیان گروه مهندسی آب. دانشگاه اراک. 

(۲)نویسنده‌ی مسئول: استادیار. گروه مهندسی آب. دانشگاه بین‌المللی امام خمینی(ره). 
(۳) استادیان گروه مهندسی آب: دانشگاه اراک. 


(۶) کارشناس: گروه مهندسی آب دانشگاه اراک. 


۷۸ 


مقد مه 

اکثر رودخانه‌ها به‌صورت پیچانرود هستند و دارای 
جریان ثانویه‌ی قوی می‌باشند که سبب رسوبگذاری در 
ساحل داخلی و فر سایش در ساحل خارجی می‌گردد. 
از طرفی دهانه‌ی آبگیر در پیچ رودخانه بایستی در 
مکانی قرار بگیرد که حداکثر میزان قدرت جریان انویه 
در آنجا باشد. این مسأله سبب می‌گردد که رسوبات در 
باتهم قاری بش #رباتهل دا کی سر کته کشزدق 
کمترین میزان رسوب وارد د ها نه‌ی آبگیر شود. 
هم‌چنین» تغییرات در پیچ رودخانه‌ها و فرسایش دیواره 
و بستر آنهاه سبب حرکت و مهاجرت پیچ رودخانه و 
تخریب تأسیسات اطراف آن و اراضی کشاورزی و 
ایستگاه‌های پمپاژ مجاور می‌گردد که ذشان‌دهنده‌ی نیاز 
به تایه ال کر مر وان و مکان سش ته تیان تاکز نهر 
تنش برشی در پیج به‌وسیله‌ی مدل‌های ریاضی است. 

شتار )٩۳6/121(‏ و مورتی (۷10:05) [1] در سال 
۲ برای شبیه‌سازی عددی جریان در حالت آشفته 
از معادلات رینولدز متوسط گیری شده در عمق برای 
حالت دوبعدی استفاده کردند. صفرزاده گندشمین 
و صالحی نی شابوری [۲] در سال ۱۳۸۵ برای برر سی 
الگوی جریان در آبگیر جانبی از نرم‌افزار فلوث" ۰ 
استفاده نمودند. ابشان برای شبیه‌سازی آشفتگی از 
مدل‌های دو معادله‌ای ع- استاندارد و 1-0 استفاده کرده 


و عملکرد مدل 0 را در پیش‌بینی نواحی گردابه‌ای و 
جدایی جریان بهتر از مدل ع-1 ارزیابی کردند. امیدبیگی 
سه‌بعدی جریان در آبگیر جانبی پرداختند. ایشان برای 
مدل سازی آش_فتگی از مدل‌های ع-ع هو 135۷ 
استفاده کردند. مقایسه‌ی نتایج تفویت آمله ۱۱ ملاع 
و با نتایج آزمایشگاهی نشان داد که الگوی جریان 
داخل کانال اصلی‌به‌وسیله ی هر سه مدل به‌خوبی 
پیش‌بینی شده اما در کانال آبگی مدل ]67۷ نسبت به 
دو مدل دیگر بهتر عمل کرده و قابلیت مدل 6-0 بیشتر 
از مدل ع- می‌باشد. شاملو و همکاران [4] در سال 


۹ به بررسی ناحیه‌ی جدایی جریان در آبگیر جانبی 


نشریه‌ی مهندسی عمرآن فردوسی 


مقایسه‌ی کاربرد مدل‌ها یآشفتگی 195 و گت در پیچ تند 


پرداختند. ایشان براساس نتایج به‌دست آمده گزارش 
نمود ند که مدل آشفتگی 15 به علت عدم اعمال 
محدودیت هموژنی لزجت گردابه‌ای» در پیش‌بینی 
الگوی جریان بهتر از سایر مدل‌های به کار رفته عمل 
می‌نماید. هان (۲10) [5] در سال ۲۰۱۰ به برر سی سه 
مدل ععل 0( کت و 88 در یک پیچ تند 
پرداخت و نتیجه گرفت که مدل 85۷ تطابق مناسبی با 
مقادیر آزمایشگاهی دارد. ون بالن (معاظ 2 ۷) 
و همکاران [6] در سال ۲۰۰۹ به برر سی مدل 1۳28 در 
پچ ملایم کا نال پرداخت ند. هم‌چنین تحقم مات بویج 
([ز30) [7]. درگاهی [8]» کانگ و چوی [9] و 
خسرونژاد و همکاران [10] نیز نشان دادند که مدل ع-1 
دارای محدودیت‌هایی در پیش‌بینی صحیح میزان جریان 
انویه است و تخمین ضعیف‌تری از ناحیه‌ی جداشدگی 
می‌نماید. از اين‌رو, به‌نظر می‌رسد برای شناخت کامل 
الگوی جریان در پیچ‌ها؛ بایستی آزمایش‌های بسیار 
زیادی را در آزمایش‌گاه و در پیچ های با اندازه‌های 
متفاوت انجام داد که انجام آن به سادگی مي سر نيٍ ست. 
بتابزاین .هدف اصتلی انجام این پژوهت» برزنستی 3 
مقایسه‌ی توانایی مدل‌های عددی آ شفتگی در شناعت 
الگوی جریان در پیچ‌های تند می‌باشد تا در صورت 
امکان با کمک آن‌ها بتوان با آزمایش‌های بسیار کمتر و 


برای شرایط مختلف. الگوی جریان را پیش‌بینی نمود. 


مواد و روش‌ها 

مطالعات آزمایش‌گاهی این پژوهش در آزمایشگاه 
هیدرولیک دانش‌گاه ,۳۳1 لوزان سوئیس, بر روی 
فلومی از جنس پلک سی گلاس, با زاویه‌ی مرکزی ۱۹۳ 
درجه انجام گردید. شعاع خط مرکزی ۶ برابر ۱/۷ متر 
می‌باشد که بیانگر نسسبت 1/3 برابر با ۱/۲ متر و 
نشان‌دهنده‌ی یک خم با پیچ تند می‌با شد. کف کانال از 
رسوباتی به قطر ۲ میلی‌متر پر شده است که پس از سه 
هفته تزریق رسوب و ایجاد توپوگرافی توسعه یافته. 
برای انجام آزمایش‌ها منجمد گردیده است. جدول (۱) 
خلاصه‌ای از شرایط آزمایشگاهی را نشان می‌دهد. 


سیالن بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


جواد مظفری- امیر صمدی- سید اسدالله محسنی موحد- داود داودمقامی 


جدول ۱ مشخصات جریان در فلوم آزمایشگاهی 


نسبت نست نسبت شعاع 
رینولدز ۱ ۲ ۲ عدد 
۱ عرض به | شعاع قوس | قوس به ۱ 
برشی ۱ ۱ ۱ فرود 
عمی به عمی عرعن 
*2 8/1 1/1 :11 39 
0 6 6 6 0 
۱۸۷ ۱۳ ۷ ۱/۳۱ ۰/۵ 


۷۹ 


: 0 بر کت عمق 
عدد شیبت 
۱ برشی متوسط | متوسط | دبی 
ریولدر برری 1 نان وم 
16 ۹۶ 0 10 1۳ 0 
(8 8 (5-1) (5-1) هس _ | (۱۶-۱) 
۰ | :۰/۰۰ 120 ۹ ۹۸ ۳ 


شکل ۱ شبکه‌ی محاسباتی ایجاد شده در نرم‌افزار گمبیت 


مدل‌های مورد استفاده در نرم افزار فلوئئت 
نرم‌افزار فلوژنت برای اجرا و شبیه‌سازی الگوی 
ستبهلی جزریان فر باقت ند [11] عدل هش سعادله‌ای 
تنش رینولدز (8[۷) از سری معادلات ناویر- 
| ستوکس با رینولدز متو سط زمانی (۵(5؟) می‌با شد 
که در این پژوهش به‌کار رفته است. مدل شبیه‌سازی 
گردابه‌ی بزرگ (1۳8) نیز با مدل‌سازی گردابه‌های 
بزرگ» الگوی جریان را تعیین می‌کند. بنابراین مدلی با 
هزینه‌های اجرایی کمتر از روش شبیه‌سازی عددی 
م ستقیم (10(5) و هزینه‌های اجرایی ب, شتر از مدل‌های 
۹ می باشد. مدل 15 یک مدل فیلتر را برای 
نشان دادن رفتار پیچه های کو چک فرض می‌کند و 
ساختار پیچه‌های بزرگ سه‌بعدی وابسته به زمان را در 
جریان آشفته محاسبه می‌کند. 

برای اجرای نرم‌افزار فلوئنت ابتدا بایستی شبکه‌ی 
جریان در کانال مورد نظر ساخته شود. برای شبکه‌سازی 
از نرم‌افزار کمبیت (02۳001)) استفاده شده است. در 
این تحقیق ٩۱۰‏ هزار گره محاسباتی در نرم‌افزار تشکیل 
ايی بع کی ساسا هاپس 6 


سال بیست و هفتم» شماره‌ی یکه ۱۳۹۶ 


چندین شبکه‌ی مختلف در نرم‌افزار گمبیت و سپس 
فلوئنت و با کوچک کردن شبکه تا حد ممکن به‌د ست 
آمد. در نزدیک دیواره» شبکه‌ی محاسباتی تا اندازه‌ای 
که در نرم‌افزار گمبیت امکان دا شته و ایجاد حجم منفی 
نکرده. ریز شده است. مقدار عبارت ۷/ م>+17 بزرگ‌تر 
از ۷۰ می‌باشد و در نتیجه جریان کاملاً زبر اسست و 
زیرلایه‌ی لزج وجود نخواهد داشت. شبکه‌ی محاسباتی 
به‌صورتی انتخاب شده بود که تا حد امکان +۷ در 
محدوده‌ی لا به‌ی میانی (1876 ۲6۲/ظ) قرار نگرفته 
باشد و بیشتر از ۳۰ شود. 

کل (۱) کانال را با فسبکه‌قی تحانسباتی, تفنان 
می‌دهد. قابل ذکر است که کف کانال از یک توپوگرافی 
توسعه‌یافته تشکیل شده است و بنابراین اندازه‌ی دیواره 
خارجی بیشتر از دیواره‌ ی داخلی در شسبکه‌ ی 
اندازه‌گیری می‌باشد. در این شبکه‌ی محاسباتی ۱۵ گره 
در عرض کانال ۳۵ گره در ارتفاع و ۰ گره در طول 
کانال ایجاد شده‌است. در گره‌های موجود در عرض 
کانال با نزدیک‌تر شسدن به دیواره‌ی جانبی. فاصله‌ی 
گره‌ها نیز کاهش یافته است. 


نشریه‌ی مهندسی عمرآن فردوسی 


برطبق روابط ون راین (مزن ۷۵۵) [12] در سال ۱۹۸۶ 
در نظر گرفته شده است. طبق رابطه‌ی ون راین» زبری 
بستر برابر است با : 

)۱( (/۵*-1.1۵0+ م34 < مک 
دیون‌ها می‌باشد. هم‌چنین ۸27.38 می‌باشد که ظ 
ترانن باعمی انب است. ارتفاع دیون‌ها نیز از رابطه‌ی 


زیر به‌دست ی 


03 
)۳( (25-7( 031 1-6) ال ِ 


وه - تم _ 


)۳( 2 1 
(+۷) 
که در آن, ,۱۳ سرعت بر شی بحرانی می‌با شد که 
7 ۷ 
از دیاگرام شلدز به‌دست می‌اید و 2 بدا 


می‌با شد. لا نیز سرعت متو سط جریان و ") ضریب 
شزی ذره می‌باشد که برابر است با: 


12۳ 
0-1102 ۶ 63 


که در اين رابطه. م6[ شعاع هیدرولیکی می‌باشد. 
شرایط مرزی استفاده شده نیز عبارتند از: 
۱- ورودی (1000): برای ورودی شرط مرزی -171255 
اعلس-1109 يا دبی جرمی ورودی در نظر گرفته شده 
است. 
۲- خروجی (0۷10۷0): شرط خروجی -۳۲695۲6 
6 قرار داده شده است. 
۳- دیواره (۷۷2[1): برای دیواره‌ی خارجی و داخلی که 
از جنس پلکسی گلاس می‌باشند. شرط مرزی ۷911 
درنظرگرفته شده و زبری آن برابر زبری ۳۷ می‌باشد. 
زبری کف نیز برابر ۰/۰۰۷ قرار داده شده‌است. 
- سطح (1۵60ا9): برای سطح آب» شرط مرزی 
77 قرار داده شده‌است. بدین معنی که سرعت 
عمودی در سطح آب. صفر درنظرگرفته می‌شسود. 
نکته‌ای که در ان خاب این شرط مرزی مورد نظر 
می‌باشد. این است که در صورتی این شرط قابل 


نشریه‌ی مهندسی عمرآن فردوسی 


مقایسه‌ی کاربرد مدل‌ها یآشفتگی 1951 و گت در پیچ تند 


استفاده است که سطح آب مشخص باشد و مانند یک 
دیواره در نظر گرفته شود. در پیچ‌ها به‌دلیل تغییرات 
موجود در سسطح آب که توسط نیروی گریز از مرکز 
ایجاد می شود سطح آب در دیواره‌ی خارجی بالاتر از 
دیواره‌ی داحلی می‌باشد و درنظر گرفتن یک سطح آب 
با ارتفاع ثابت می‌تواند خحطا ایجاد کند. به‌هر حال در 
تحقیقاتی که توسسط ون بالن و همکاران [6] بر روی 
پیچ‌های ملایم و تند انجام شد. بیان گردید که مقدار 
این خطا ناچیز می‌باشد و تأثیر آن بر روی الگوی جریان 
بین بردن تعطای مذکور در شبکهی محاسیاتی؛ از سطح 
آب اندازه گیری شده اسخفاده کر تن بدین تردیب 
شبکه‌ی محاسباتی در فضای واقعی بین توپوگرافی 
توسعه‌یافته و سطح واقعی آب طراحی شده‌است. 


نتایج و بحث 

کمیت‌های متوسط گیری‌شده در عمق معمولاً برای 
شریح سیمای عمومی میدان جریان و توزیع سرعت 
در جریان‌های سهبعدی پیچیده مورد استفاده فرار 
می‌گیرند. شکل (۲) سرعت طولی متوسط گیری‌شده در 
عمق (ل1) را که با سرعت متوسط جریان نرمالیزه شده 
است (09/17) نان می‌دهد. جریان تا حدود مقطع 
عر ضی ۳۰ درجه دارای جدا شدگی از دیواره‌ی داخلی 
نمی‌باشد. منطقه‌ی جداشدگی بین ۳۰ و ۱۲۰ درجه 
رت گنه یازع شته عرش وروی ۲0 مرس ات 
که در حدود 1۰ درصد عرض کل می‌باشد. 

شکل (۳) تغییرات سرعت عر ضی متو سط‌گیری 
و هی هی را که مسر یسوط زان 
پمال و کته اشست: (ا )نان ی دهاز سوه 
مقطع عرضی ۱۵ درجه تا مقطع عرضی ۷۵ درجه 
مقدار متوسط گیری‌شده‌ی سرعت عرضی مثبت است و 
به سمت ساحل خارجی می‌با شد. از مقطع عر ضی ۷۵ 
در جه تا نزدیک ۱۳۵ در جه عمد تا مقدار سرعت 
عرضی متوسط گیری‌شده در عمق, منفی و رو به ساحل 
داخلی می‌با شد. از مقطع ۵ درجه به سمت خروجی 
حم» دوباره لین مقدار مثیت گردیده است. 


میالن بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


جواد مظفر ی - امیر صمدی- سید اسداله محسنی موحد- داود داودمقامی ۸۱ 
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شکل ۲ تغییرات سرعت عرضی متوسط گیری‌شده در عمق (170/1)- براساس مدل آزمایشگاهی 


وی._.-. 


شکل ۶ سرعت طولی متوسط گیری شده در عمق (ا/1[5)- پیش‌بینی شده با مدل 115۷ 


بررسی و مقایسه‌ی الگوهای متوسط گیری شده از نرم‌افزار فلوئنت به‌دسست نمی‌آیند. برای تعیین این 
الگوهای متو سط گیری شده در عمق به صورت هستقیم الگوها در ابتدا پروفیل‌های سرعت طولی در 


سال بیست و هفتم. شماره‌ی یک ۱۳۹۶ نشریه‌ی مهندسی عمران فردوسی 


۸۲ 


موقعیت‌هایی از کانال که در آزمایشگاه نیز اندازه‌گیری 
شده بود به‌دست آمد. این موقعیت‌ها عبارتند از ۱۵ 
مقطع که در هر مقطع عرضی ۰ پروفیل سرعت 
استفاده شده‌ا ست. در مرحله‌ی بعد. از تقاط موجود در 
هر پروفیل سرعت میانگین گیری وزنی شد تا مقدار 
مقوشیط گیرع هی ان فرعم بادست ابید تن اب 
مقادیر وارد برنامه‌ی متلب (۷120120) شد و الگوهای 
متوسط‌گیری شده در عمق رسم شدند. شکل (1) 
سرعت طولی متو سط گیری شده در عمق را که تو سط 
مدل هفت معادله‌ای تنش رینولدز پیش‌بینی شده‌است 

با توجه به شکل (4) هم‌چنان ناحیه‌ی بیشینه 
سرعت در مقطع بین ۰ تا ۳۰ درجه متمایل به‌سمت 
داخلی قوس پیش‌بینی شده است. اما همان‌طور که قبلا 
بیان شد ناحیه‌ی بيٍ شینه سرعت اندازه‌گیری شده در 
مرکز قرار دارد. بین مقاطع عرضی ۳۰ و ۷ در جه 
حالی که مقدار اندازه گیری شده به‌سمت ساحل خارجی 
نجل تست گر هم‌چ نان در مقطع عرضی ۰ در حه 


پیش‌بینی شده‌است که تا حدودی به سمت زاویه‌ی 


مفایسه‌ی کاربرد مدل‌ها یآشفتگی 1 و ٩‏ در پیچ تند 


قوس ۸۰ در جه نیز حر کت کرده‌است اما هنوز با 
منطقه‌ی جداشدگی اندازه‌گیری شده کاملاً متفاوت 
است. هم‌چنین منطقه‌ی جدا شدگی اندازه‌گیری شده تا 
خحط مرکزی کانال پیشروی کرده‌است اما مقدار 
پیش‌بینی شده در مدل تنها ناحیه‌ی کوچکی را در بر 
نی گیرکد از مقطع عرضی ۰ تا ۱۳۰ در جه از حط 
مرکزی تاساحل خارجی. پیش بینی و اندازه‌گیری 
یکسان می‌باشند اما از خط مرکزی به سمت ساحل 
داخلی مقدار پیش‌بینی شده همان ناحیه‌ی جداشدگی را 
شتان می‌د هد که متفاوت با اندازه گیری می‌باشد. از 
مقطع عرضی ۰ درجه تا مقطع عرضی ۰ درحه. 
مقدار بر شینه سرعت پیش‌بینی شده به سمت ساحل 
خارجی می‌باشد و تا خروجی پیچ در ساحل خارجی 
جای گرفته‌است اما مقدار اندازه‌گیری شده‌ی آن تا 
حدودی متفاوت است و از مقطع عرضی ۰ درحجه در 
ساحل خارجی جای می‌گیرد. 

شکل (۵). سرعت طولی متوسط گیری‌شده در 
عمق را که با مدل 8 پیش‌بینی شده‌است نشان 
می‌دهد. در این شکل نیز در مقطع ۰ تا ۲۰ درحه هسته 
ببرشینه سرعت در سمت ساحل داخلی قرار دارد اما تأ 
حدودی به سمت مرکز نیز متمایل شده است. 


شکل ۵ سرعت طولی متوسط گیری‌شده در عمق (1[9/0)- پیش‌بینی‌شده با مدل 1۳5 


نشریه‌ی مهندسی عمرآن فردوسی 


سا بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


جواد مظفری- امیر صمدی- سید اسدالله محسنئی موحد- داود داودمقامی 


بین مقاطع عرضی ۳۰ تا 1۰ در جه. پیش بینی 
جریان توسط مدل 8 نشان می‌دهد که هسته‌ی 
مرکزی بیشینه سرعت به سمت ساحل خارجی حرکت 
کرفباابیت هیا معاکی اندانه گر فده قطان مت تیش 
را نشان می‌دهد. بیشینه‌ی منطقه‌ی جداشدگی جریان که 
توسط مدل پیش‌بینی شده‌است تقریباً در مقطع ۷۵ 
درجه قرار دارد که نزدیک به مقدار پیش‌بینی شده‌است. 
این در حالی است که بیشینه‌ی منطقه‌ی جدا شدگی در 
متا ۱۹4 کین اش عیشت نم ار 
شده‌است در مقطع عرضی ۸۵ درجه قرار گرفته‌است و 
با مقادیر اندازه‌گیری‌شده تطابق ندارند. اما منطقه‌ی 
جداشدگی در مدل 1۳7٩‏ متمایل به سمت مقطع 
غرضسی ۷۹ دز جه شسده‌اسست و تفریبا با مفادیز 
اندازه‌گیری‌شده تطابق پیدا کرده‌اسست. از مقطع ۱۵۰ 
درجه به سمت خروجی پیش‌بینی مدل نشان‌دهنده‌ی 
قرارگیری هسته‌ی بیشینه در ساحل خارجی می‌باشد که 
انش قس 6 انتا رها ای نمی زان 
می‌د هد. از مقطع ٩۰‏ تا ۱۲۰ در جه. مقادیر پیش‌بینی 
شده و اندازه‌گیری‌شده تقریباً یکسان می‌باشند و 
هسته‌ی بيشینه سرعت در هر دو الگوی جریان بین خحط 
مرکزی و ساحل خارجی قرار دارد. در بین مقاطع ۱۲۰ 
تا ۱۵۰ درجه نیز الگوی پیش‌بینی‌شده و اندازه‌گیری‌شده 
تطابق مناسبی را نشان می‌دهند. هسته‌ی بیشینه‌ی جریان 
هنوز در نزدیکی خط مرکزی و متمایل به ساحل 
خارجی است که با نزدیک شدن به مقطع ۱۵۰ درجه تا 
حدودی به ساحل خارجی نزدیک‌تر می‌گردد. 

مقایسه‌ی سرعت طولی متوسط گیری‌شده در عمق 
که به‌و سیله‌ی این دو مدل پیش‌بینی شده است ذشان 
می‌دهد که در ناحیه‌ی جداشدگی جریان» مدل 115[۷ 
جریان برگشستی آب را پیش‌بینی نکرده‌است و حداقل 
تا ان وفع اما هل سره رب 
را ذشان می‌دهد که ذشان‌دهنده‌ی بو سس آب و وجود 


گردابه در اين منطقه می‌باشد. البته ناحیه‌ی مذکور 


شا بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


۸۳ 


به‌دلیل عمق کم در آزمای شگاه اندازه‌گیری ذشده است. 
شکل (0) نیز سرعت عر ضی متو سط‌گیری شده در 
عمق را که توسط مدل تنش رینولدز پیش‌بینی شده 
است نشان می‌دهد. 

مقدار بیشینه‌ی سرعت عرضی متوسط گیری شده 
در عمق در مقطع ۰ درحه اتفاق افتاده ا ست که مقدار 
بیشتر از مقدار اندازه‌گیری شده 
می باشد ول هر دو م قدار در یک موقعیت اتفاق 
افتاده‌ا ند. از مقطع ورودی تا حدود مقطع عرضی ۷۵ 
درجه مقادیر اندازه‌گیری‌شده نشان‌دهنده‌ی سرعت 


عر ضی متو سط‌گیری شده‌ی مثبت و به سمت ساحل 
خارجی می‌باشد و پس از آن تا مقطع عرضی ۱۳۵ 
درجه مقدار آن منفی و به سسمت ساحل خارجی 
می‌با شد و پس از آن دوباره به سمت ساحل خارجی 
تغییر جهت می‌دهد. در مقادیر پیش‌بینی‌شده توسسط 
مدل. سرعت عرضی متوسط گیری‌شده در عمق از 
ورودی تا حدود مقطع عرضی ۸۰ درجه به سسمت 
ساحل خارجی می‌باشد و سپس مقدار آن منفی شده و 
تا حدود مقطع عرضصی ۰ درجه به سمت ساحل 
داخلی می‌باشد. پس از آن دوباره به مسمت ساحل 
خارجی برمی‌گردد و مقدار آن مثبت می شود. مقایسه‌ی 
سرعت عر ضی متو سط‌گیری شده در عمق بین مدل 
مذکور و مقدار اندازه‌گیری شده ذ شان‌دهنده‌یی اختلاف 
تقریتاً ۵ درجه‌ای می‌باشد. در مقاطع عرضی ٩۰‏ تا ۱۲۰ 
درجه هسته بیشینه سرعت عرضصی در الگوی 
«اشاهله‌ای نزدیک یه مانعل خاویخی | سظه درخعالی که 
در مدل محاسباتی نزدیک به ساحل داخلی می‌باشد. لذا 
این مقادیر صسحیح پیش‌بینی نشده‌اند. مدل مذکور 
الگوی جریان را به‌صورت کلی صسحیح پیش‌بینی 
کرده‌ا ست اما در پیش‌بینی مقادیر بٍ شینه و موقعیت آن 
دارای عدم تطابق هایی با مقادیر | ندازه گیری‌شده 


ی 


نشریه‌ی مهندسی عمرآن فردوسی 
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.--س۳ 


شکل 71 سرعت عرضی متوسط گیری‌شده در عمق ([110/1)- پیش‌بینی‌شده با مدل 115۷ 


شکل ۷ سرعت عرضی متوسط گیری شده در عمق ([11/1)- پیش‌بینی‌شده با مدل 125 


شکل (۷) سرعت عرضی متو سط گیری شده 
در عمق را که با مدل ٩‏ پیش‌بینی شده‌است نشان 
می‌دهد. این مدل همانند مدل پیشین دارای پیش‌بینی 
تقریبا مشابهی برای الگوی سرعت عرضی 
متوسط گیری‌شده در عمق می‌باشد. در بین فاصله‌ی 
مقاطع ۰ تا ۱۲۰ درجه. دو سرعت عر ضی بٍ شینه را 
تفن رسک کزان رفک مان انار 
اندازه گیری شده می‌با شند اما مقدار بیشینه‌ی نزدیک به 
ساحل داحلی بیشتر از مقدار اندازه گیری شده پیش‌بینی 
شده‌است. 
الگوی پیش‌بینی‌شده‌ی سرعت عرضی از مقطع 
ورودی تا حدود مقطع عرضی ۵ در حه مابت و 


نشربه‌ی مهندسی عمرآن فردوسی 


به سمت ساحل خارجی می‌با شد و این روند تا حدود 
مقطع ۱۳۵ در جه ادامه دارد. پس از آن الگوی جر یان 
منفی شده و به سمت ساحل داخلی قرار گرفته ااست. 
این پیش‌بینی روند بهتری را نسبت به مدل پیشین نشان 
می‌د هد. در مدل قبل اختلاف تقریباً ۵ در جه‌ای بین 
الگوی | ندازه گیری شده و پیش‌بینی شده وجود 
داشت که اپن اختلاف در مدل 1۳28 از بین رفته است. 


مقایسه‌ی پروفیل‌های سرعت طولی اندازه گیری شده 
و پیش‌بینی شده 
انتخاب شده م شخص کرده ا ست. در این حدول ابتدا 


سالن بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


جواد مظفری- امیر صمدی- سید اسدالله محسنئی موحد- داود داودمقامی 


مقدار خطای چهار پروفیل سرعت طولی انتخاب‌شده از 
هر مقطع عرضی. تعیین گردیده و سپس متوسط خطای 
آنها متشی شاواستی بطرر رشان فضای مین 1۳ 
در مقطع عرضی صفر درجه و موقعیت 5 ۱/۲ 
در صد می‌با شد. هم‌چنین متو سط خطای سرعت طولی 
در مقطع ورودی برای پیش‌بینی مدل نان 1/۵ در صد 
است که کمتر از خطای مدل 11801 می باشد. 8 
نشان‌دهنده‌ی جهت شعاعی می‌باشد. به‌طور مثال 21.35 
به معنای فاصله‌ی ۱/۳۵ متری از مرکز پیچ می‌باشد و با 
توجه به اين که دیواره‌ی داخلی با مرکز پیچ ۱/۰۵ متر 
فا صله دارد. 01.35 ذشان‌دهنده‌ی موقعیت ۰/۳ متری از 
دیوار داخلی می‌باشد. هم‌چنین 11.7 نشان‌دهنده‌ی 
فاصله ۱/۷ متری از مرکز پیچ می‌باشد و موقعیت آن 
سط پر کی کانان اس 09 قوش تیش ۷سا ان 
دیواره‌ی خارجی و 92.2 موقعیت ۰/۱۵ متری از 
خر ار ار شا ای ات 


حدول ۲ درصد اختلاف بین مقادیر پیش‌بینی شده و اندازه‌گیری 


شده پروفیل سرعت طولی 
زاویه 
مرکزی | مدل | ک3له | 17« | ک0جه | 22« | میانگین 
(درجد) 
11۳5 ۸ ۸/۲ ۳ | ۰/۶ 1/۵ 
۱ | ۰/۶ ۱ ۱۱/۲ ۱ ۱۱/۶ | 1/۲ 1 ۰ ۷/۳ 
وق | ۱۲۸ | ۵ | ۵ | ۶۸ 1۵ 
ِ | ۱۲۸۲ ۱ ۱۲ ] ۹/۸ ۱ ۵/۸ ۱۰ 
و | از ۱۳ ال ۱ ]۹ 1 ۰ ۷۸۷ 
| ۱۳ ۱۱ ۱۷ ۱۳/۷ 
و | | ۹۸۷ | ۸۸ 1 ۱۵ 1 ۱۲/۵ 
0 | ات لا با ۲۵ | ۹ ۱ ۷۸۵ 
ویر | ۸۸۶ ۱ ۱۹ ]| ۱۸۷ | 1۱۸ ۷۸۷۰ 
ت ۸ | ۲۰ 1 ۰ | ۶/۷ | ۱۰/۲ 
ور ۱ ۱۲ ۸ ۳ | 1/۸ ٩/۵‏ 
۳ 7 ۱۰ ۹ ۲ ۸/۲ ۱ ۱۳/۸ 
و | ۶/۶ | ۱۱/۵ ۵ 2/۵ 2/۹ 
ِ ۸ | ۲۵ | ۷/۵ | ۵۸۷ | ۳۸۷ | ۱۰/۵ 


بر طبق جدول (۲) مدل 125 با ۷/۹ درصد خطای 


میانگین کل مقاطع عرضی مقایسه شده بهترین نتایج را 


سال بیست و هفتم. شماره‌ی یک ۱۳۹۶ 


۸۵ 


دربرداشته است. مدل 491 نیز دارای ۱۰/۶ درصد 
خطا می‌با شد. در هفت مقطع عر ضی انتخاب شده تنها 
در زاویه‌ی مرکزی ٩۰‏ درجه مدل 1۳5 خطای بیشتری 
را داراست؛ البته باید به این نکته توجه داشت. تنها 
خطای پیش‌بینی ۲۱ موقعیت مکانی در کانال بررسی 
شده است و اگر در زاو یه‌ی مرکزی ٩۰‏ در جه 
موقعیت‌های دیگری انتخاب می‌شدند امکان تغییر نتایج 
وجود داشت و احتمالاً با افزايش تعداد پروفیل های 
سرعت طولی در اين مقطع. مدل 125 جواب بهتری 


می‌داد. 


نتیجه گیری 
برای بررسی مدل های عددی آشفتگی جریان و 
مقایسه‌ی نتایج آنها با داده‌های واقعیء از داده‌های یک 
آبراهه با پیج تند و توپوگرافی توسعه‌یافته اسستفاده 
گردید. مطالعات آزمایشگاهی بر روی فلومی با شعاع 
انحنای ۱/۷ متر و با زاویه‌ی مرکزی ۱۹۳ درجه انجام 
شده است. جریان ۱۳ لیتر بر ثانیه بر روی ر سوبات با 
قطر ۲ میلی‌متر اجرا شد که پس از این که توپوگرافی 
کاملاً تو سعه یافت. برای انجام اندازه‌گیری‌ها ر سوبات 
منجمد شد. برای تعیین الگوی جریان به‌وسیله‌ی 
مدل‌های آ شفتگی از نرم‌افزار فلوئنت و گمبیت استفاده 
گردید و مدل‌های آشفتگی 153۷ و ۳٩‏ مورد برر سی 
قرار گرفتند. بررسی‌ها نشان داد که هر دو مدل تا 
حدودی الگوی کلی جر یان را پیش‌بینی کردها ند اما 
بخش‌های مهم جریان به‌وسیله‌ی مدل 1551۷ به حوبی 
پیش‌بینی نشده بود. مدل 181۷ نتوانسته بود منطقه‌ی 
جدا شدگی جریان را به صورت مناسبی پیش‌بینی کند و 
به‌جای موقعیت زاویه‌ی ۷۵ درجه زاویه‌ی ۸۵ درجه را 
بعنوان متطقه‌ی: خداشتد کین پیش بیش کراده‌اسستتا: ما 
مدل 15 موقعیت این ناحیه را به‌صورت مناسسبی 
پش‌بیتی. گر همین هن پیش بیش موته بت مسبت که 
پیشینه‌ی سرعت طولی: مناسب‌ترین پیش‌بینی برای 


مدل 25 بوده است. در سرعت عرضی متو سط گیری 
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۸ 


شده در عمق نیز مدل 178 تغییر جهت جریان عرضی 
و موقعیت آن را بهتر پنش‌بینی کرده است: نهایتً تعطای 
میانگین میزان اختلاف بین پروفیل‌های پیش‌بینی شده و 
اتداژه گر له بر از سل فا کر از ۸ درض تا 
به‌دست آمد و این مدل دارای کمترین میزان خطا در 


پیش ‌بینی بوده‌است. 


ات 


مقایسه‌ی کاربرد مدل‌ها یآشفتگی 1851 و گت در پیچ تند 


تشکر و قدردانی 
این پژوهش از محل طرح پژوهشی داخلی دانشسگاه 
اراک با موضوع «بررسی مدل‌های آشفتگی در پیچ 
رودخانه با کمک نرم افزار ۳11/11 تامین اعتبار شده 
تن پوریستتند کان بر خود لازم می‌دانند که از معاونت 


پژوهشی و فناوری دانشگاه اراک قدردانی نمایند. 


9۲۳۵ولر و "واعصصعط صووم صا م1 ۵۶ وونل ۵۶ بای آهم1تمصصتاه ۵ ویکلیک روطاسا مه رفظ رتمتامطه ...1 
(1996) ,651-675 .00 و5 ٩0,‏ م34 ۷۵۱۰ رع‌صه۱عظ۱۵0 مط 1 هب۲0۵۵۱۲۵5۵0۲ 


۲ صفرزاده گند شمین, ا. و صالحی نیشابوری ع... "مطالعه‌ی عددی الگوی جریان آشفته و بررسی کیفی انتقال ر سوب و فرسایش در 


آبگیری جانبی از رودخانه ‏ مجله‌ی فنی و مهندسی مدرس. شماره‌ی ۰۲۵ صص. ۰۱-۱۸ (۱۳۸۵). 


۳ امید بیگی. م.ع "مطالعه‌ی آزمایشگاهی و شبیه‌سازی عددی سه‌بعدی الگوی جریان در آبگیری جانبی از رودخانه در حضور صفحات 


مستغرق " پایان‌نامه‌ی کارشناسی ارشد. دانشکده‌ی کشاورزی, دانشگاه تربیت مدرس. (۱۳۸۸). 


ع شاملوی ح.» و پیرزاده» ب» "بررسی اثرات پارامترهای هندسی و هیدرولیکی بر ابعاد ناحیه‌ی جدایی جریان در آبگیرهای جانبی به‌ کمک 


نرم‌افزار فلوئنت " نشریه‌ی مهندسی عمران و نقشه‌برداری» صص. 6۹۵-1۸۷ (۱۳۸۹). 
,۱10108 ۵۶ ,1۱۳۵۵۲ .ع1وعظ]. .ملظ رتعصوها امصصقطه هون 90 فصامه ماه گم ومتاوزهامهتهطن ریقنو رصع .5 


.(2010) و۵106060) باهه‌صم رنه نصا هت4تممصمن رمصتهم م۴ وممصصمز ص۴۳ قمع 0 


0۵ نان اتمه ۵۶۴ صمتاهاناصزه و600 معتصا ریک راتمهیامصهاظ مه ر.ل ۷۷۰۹ راهه هار۱ و۷۷۰ رصعلوظ مه۷ .6 
,(2009) ,413-442 .00 ,630 ,۷۵1 ,,/۱۷۲6۵ ت۲۱ .ل ,1100۲ امصصعهم 


بل و کعصتا ۵0تته موی صز قوماگ قط گم فصممتاهانصته برع معععا مه کامممم‌تنافهع]۷( رب رژا00ظ .7 


2003(۰) 6-16 .00 با .0 و4 ۷۵۱ رع6۳6 71۱۳ 


هط رصهب هلک عز۱ مط ا اعممورقصهه کمم‌حصطرزلعو 4ج مصتاام0مصه ما لمممتفممصوتل-عمط 1 ربظ رتطهع‌تهنا .8 


2004(۰) ,852 -821 .00 ,7 ۱0۰ و29 ۷۵۱۰ وکن0]0صقا قمع وعووعع۲۵ظ ومه]ا8 


۰ .[ ,۲۳۵۱ ,۶10۷7۲ امصصهطمصوون تقانم‌ممامع۲ ۵۴ مصاام0مصه وومتاو 05مصوم۴" ریق رتمطل 220 وا رعطقک ,9 
.(2006) ,1319-1334 .00 ب11 .0 و51 ,۷۵۱ رد۳۱ 1۷۲6۰ 
عصنام0 مصط آممزمصناطر مرو ولیک و۵۳۳[ مه وق شمش رتتاممامطونوها ر.نلیت رقتصصفگ ریم ب24ز7080ومطک1 ,10 


٩0. 10, ۰‏ و133 ,۷۵۱ ر.ق عزابم۱۵۲ .7 و"عقجهها امصصعطه امیجتمطصاع1 طا کتممرعصنه کصمحصتعی صرح اما ۵۶۴ 


1123-1134, )2007( 


2006(۰) رهظ ۵۴ وماها و ۲۲1660 رتتطوم‌صها وا رمصع‌جامن ر"ع0تن وتمونا 6.3 ۳۸۲1 نابا۳۲" ور.عصا کصمیاا ۳ :11 
٩۲۵. 10 ۰‏ و110 ,۷۵۱ .۳9 0۲۰ .ل موجه 10240 ۵۵0 :1 تور راتمجقصعا کصمصتلعو وا مطزنک ۷2 .12 


1431-1456, 0984۰ 


,"1مصصعدطن ۵۶ عصوها طا ماگ اهامای ۵۶ مملاها نی آممزمحرباطر نو ریک رتصصم‌طعه۲۱ مه و۷ رلته‌2ه۱0۵ ,13 


2012(۰) ولوانا رو۱۱۵۵۲ و ۱۱6۵۱۵۲6۵ ۵۴ 6۵۵/6۲۵۴۵۵۵ ۳۱۵۵۲۳۱۵۵۵۴۵] 
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سیالن بیست و هفتم. شماره‌ی یک. ۱۳۹۶ 


